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Abstract 

We compare the gauge-fixing approach proposed by Jona-Lasinio and Parrinello, and 
by Zwanziger (JPLZ) with the standard Fadeev-Popov procedure, and demonstrate 
perturbative equahty of gauge-invariant quantities, up to irrelevant terms induced by 
the cutoff. We also show how a set of local, renormalizable Feynman rules can be 
constructed for the JPLZ procedure. 



^ e-mail: bock@physik.uni-siegen.de 
^ e-mail: maarten@aapje.wustl.edu 
^ e-mail: mco@morgan.wustl.edu 
^ e-mail: shamir@post.tau.ac.il 



1. Gauge fixing in Yang-Mills theories is well understood in perturbation theory. 
But outside of perturbation theory the situation is more complicated, because of the 
existence of Gribov copies These complications arise as soon as the theory is 
nonlinear, which is the case for nonabelian gauge theories, but also, for instance, with 
a nonlinear gauge-fixing condition in an abelian theory. 

This does not mean that Yang-Mills theories do not exist outside perturbation 
theory. It is well known that they can be defined on the lattice, using the compactness 
of the group to dispense with gauge fixing altogether. However, nonperturbative 
gauge fixing is interesting in a variety of contexts, ranging from the need to make 
contact between lattice and continuum calculations (see for example Ref. 0) to the 
definition of chiral gauge theories on the lattice (see for example Ref. [Q). For a 
discussion of the numerical implementation of the gauge-fixing method discussed in 
this paper, see Ref. 

One can try to extend the standard BRST construction to the lattice. When 
one does this, one finds, quite generally and rigorously, that the BRST gauge-fixed 
partition function vanishes [Q. This result can be heuristically explained as the result 
of pairwise cancellations of lattice Gribov copies, which occur with opposite signs of 
the Faddeev-Popov determinant. (For some recent work on overcoming this problem, 
see Refs. 0, ||.) 

A different nonperturbative method for gauge fixing has been proposed some 
time ago by Jona-Lasinio and Parrinello and by Zwanziger (JLPZ) 0, |10| (see also 
Ref. [|11]). Starting from the euclidean ungauged partition function 



Z = j PHAexp(-^i,v(A)) , (1) 

one inserts one in the form 

(2) 



_ /Pn/^exp (-^ni(A") 
^~ /PH/iexp(-5ni(A'^)) 



into the functional integral. Here and D-nh denote the invariant measures over 
the gauge field A and the group-valued scalar field /i, respectively; is the gauge 
transform of A under a gauge transformation h. S^i{A) is any gauge noninvariant 
local functional of A. Because of the invariance of I'h^ and S'inv(^), we can perform a 
gauge transformation in the numerator of Eq. (P with Eq. @) inserted, and, dropping 
a trivial factor / "Dh^, we obtain 

^ . exp(-ginv(A) -gni(^)) , . 

" SVnhe^Y>{-Sni{A^^)) ' ^ ^ 

This procedure is completely rigorous on the lattice, and sources coupled to gauge- 
invariant operators can be added without changing the argument. We also note that 
the Boltzmann weights for both integrations contained in Eq. (^ are positive. 

While this method is conceptually very simple, much less is known about how 
this nonperturbative method works in perturbation theory, as compared with the 
standard gauge fixing based on BRST invariance. Note, however, that there are 
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some similarities between the two procedures. A noninvariant term S^i{A) is added 
to the action, and (for a suitable choice) this will make the quadratic part of the 
action invertible. The contributions of different orbits are weighted properly because 
of the group integral in the denominator of Eq. which plays a role similar to 
that of the Faddeev-Popov (FP) determinant in the usual case. In the FP case, 
however, the determinant can be expressed as an integral over ghosts, and the gauge- 
fixed action including the ghost terms is local, whereas here this is not the case. 
This is an important difference, because locality is a key ingredient in power-counting 
arguments, and thus at the heart of the usual perturbative analysis of renormalization. 

It is therefore of interest to find out whether perturbation theory can be system- 
atically developed for a JLPZ gauge-fixed Yang-Mills theory, and its relation to the 
usual FP results. This question has been addressed previously by Fachin |12|, who 
analyzed the vacuum polarization for the choice S^i{A) = tr(M^y4^) at one loop. He 
concluded that for M — > oo the transverse part of the vacuum polarization agrees 
with that obtained using the FP method; the longitudinal part vanishes for M — >■ oo, 
as in Landau gauge. This equivalence for M — >■ oo at fixed cutoff was already derived 
formally in Ref. P|, and another formal discussion appeared recently in Ref. [|13 



Here, we are interested in considering the situation where M is chosen to be of the 
same order of magnitude as the cutoff. The transverse part of the one-loop vacuum 
polarization with JLPZ gauge fixing, for example, was found to contain terms pro- 
portional to p'^Sfj^i, — p^Pu times p^ /M"^ or p'^/M'^ [|12|, so, if we choose M ~ A (with 
A the cutoff), such terms are of order 1/A^, and vanish when we take A — ^ oo. We 
present below a general argument that, if we choose M to be fixed in units of the 
cutoff, perturbation theory for a JLPZ gauge-fixed theory is equivalent to that of the 
same theory gauge-fixed using the standard FP method. We discuss in some detail 
what "equivalent" means in this context. We also show that our arguments for the 
equivalence of JLPZ and FP gauge fixing can be used to construct a set of local 
Feynman rules for the JLPZ gauge-fixed gauge theory. 

2. We begin with the JLPZ gauge-fixed partition function in the presence of a 
source for the gauge field (working in euclidean space-time): 

7. n , exp (-ginv(A) - SnijA) + J^A^) 

where Si^^{A) is the gauge- invariant classical action, and Sj^i{A) is not invariant. As 
noted before, A^ is the gauge transform of A under a finite gauge transformation. 



We will take 
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SUA) = -tT{Fp, (6) 
5ni(A) = tr(MMJ), 

where we do not indicate the integration / d^x explicitly, and M a parameter with 
the dimension of a mass. We will assume that M is proportional to the cutoff for the 
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case of regulators such as the lattice, in which case M is proportional to the inverse 
lattice spacing. However, it is sufficient for our arguments that M is chosen large 
compared to all physical scales. This assumption applies in the case of dimensional 
regularization. Also, = A^jT"-, with T"' the generators of the gauge group, with 

tr(T"T^) = (7) 
[T'',T^] = if^'^T'' . (8) 

If we now consider perturbation theory, and use a regulator without power-like diver- 
gences such as dimensional regularization, we may extend the range of integration of 
the variables from — oo to oo. The invariant measure is just the flat measure 

VnA = VA=\{ dA^ix) , (9) 

and similarly for (after parametrizing h as in Eq. (|1^) below). 

We wish to show that this gauge-fixed partition function is equivalent (in a way to 
be discussed in more detail in the next section), under certain assumptions which are 
valid in perturbation theory, to the standard FP gauge-fixed form for the covariant 
gauge d^A^ = 0. The derivation is based on two ingredients. The ffist ingredient is 
that we insert a constant into the partition function, written in the form 

constant = det(n) J Vt] 6{d^A^ - -^^v) , (10) 

where Vrj denotes the (flat) measure for t], a new field which takes values in the Lie 
algebra of the gauge group. We may choose boundary conditions such as to avoid 
the trivial zero mode of □. The second ingredient is to make use of the fact that 
the physics at scales below M is not altered by adding or changing terms which are 
irrelevant in the sense of Wilson's renormalization group. 
First, let us expand Sni{A^) in g/M, writing 

h = exp{^j^e). (11) 



We find 



where 



S^A') = ^M\A;)' - M9^d,A; + ^d,9-{D,{AW + ^ (^) ' (^2) 



{D^{A)er = d^e'' + gr'"'Ale'^. (13) 

Note that —d^D^{A) is just the FP operator for the covariant gauge. 
We then shift the vector field 

A, - A', = A,- ^d,v . (14) 
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This gives 



SAA'^) = ^M\A';Y + MA';d,^^ + -{d,vr (is) 



and 



SUA) + S^A) = S,,M') + lM'{Ay + MA'ldX + l{dXy + [j^) • (16) 

The terms |M^(A'|^)^ + MA'"'^d^rj°' + \{d^rj'^)'^ cancel between the numerator and 
denominator of the integrand in Eq. (^. The terms M6"'d^A"^ and MA"^d^r]"' (which 
is equal to —Mrj°-d^A'^^ by partial integration) may be dropped because of the 5- 
function, 5{d^A' ^), in Eq. ([T0|). 

All the 0{g/M) terms are irrelevant, and may therefore be omitted without 
changing the (renormalized) theory. Doing this, we can perform the /i-integral in the 
denominator of the integrand in Eq. (^, obtaining 

Z{J) = det{n) JV7]VA' 6{d^A'^)det'/\~d^D^{A')) (17) 

xexp (^-S-UA') + J^.A^ - ^7]n{-d^D^{A'))-'nr]^ . 

Using the fact that physical quantities do not change when we replace the source 
term J^A^ by a new source term J'^A' ^ coupling to A' ^ instead of to A^ {A^ and A' ^ 
are equivalent interpolating fields), we can now also perform the rj integral. Dropping 
the primes, we obtain 

Z{J) = jvA 6{d^A^) det{-d^D^{A)) exp {-S^A) + ■ (18) 

Representing the 5-function as 

5{d,Al) oc hm exp (^-1(9,^^^^ , (19) 

and introducing an algebraic field B and ghost fields c and c, this can be recast as 

Z{J) = Mm J VAVBVcDc (20) 

1 



xexp ( -SUA) - -^B' + tBd^A^ - c{-d^D^{A))c + J^A 



This is the standard BRST-invariant form of the FP gauge-fixed partition function 
in Landau gauge. Of course, correlation functions of gauge-invariant operators are 
independent of ^. We observe that, if we would take the limit M — > cxd before removing 
the cutoff, our argument constitutes an alternative derivation of the equivalence of 
Eq. (^) to Landau gauge given in Ref. Note that the standard way of gauge fixing 
employed in lattice QCD computations is formally equivalent to this limit. What is 
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new here, is that we do not take the hmit M — > oo first, but keep it at the order of 
the cutoff. Nevertheless, the parameter M has disappeared from Eq. (PO]). 

The derivation given above is vahd only in perturbation theory. We assumed 
that the 9- and ?7-integrals converge, i.e. that the FP operator —d^D^{A) has only 
positive eigenvalues. This is not in general the case, but it is true in perturbation 
theory. 

If we use a regulator with a hard cutoff, such as the lattice, additional subtractions 
will be needed in order to remove power-like divergences, which may appear as a 
consequence of dropping irrelevant terms. Also note that, even though Eq. is 
linear, in general the invariant measure written in terms of the Lie-algebra valued 
fields is nonlinear for such a regulator, and it would therefore change under this 
transformation. However, this nonlinearity is proportional to the coupling constant 
g, and therefore the effects of this shift are of order g/M. 

3. In this section, we will discuss in more detail what we mean by "equivalent." 
It is clear that, in general, correlation functions of the form {A'Ji^{x)A'l{y) . . .) are not 
the same in the JLPZ and FP versions. A trivial example is the 0{g^) two-point 
function {Af,{p)A,y{q)) = S{p + q)Gf,u{p), with G^,y{p) equal to 

G'^'Hp) = (21) 

in the theories defined by Eqs. (^) and (PU[), respectively. This is not in contradiction 
with the general argument presented above, because of the change of interpolating 
field, A^ A'^. 

The equality of the transverse part of the two-point function at tree level (Eq. (|21|)) 
suggests in what sense the two theories defined by Eqs. (|) and (^0|) are equivalent. 
These two theories are not identical, because of the fact that we dropped irrele- 
vant terms in going from the JLPZ to the FP version. The correct statement is 
that physical quantities (which are necessarily extracted from correlation functions of 
gauge-invariant operators) in the JLPZ version can be mapped into those of the FP 
version by a finite renormalization of the bare coupling constant g. In the FP version 
of the theory the fact that only a coupling constant renormalization is needed follows 
from BRST invariance. In the JLPZ version, the same follows from the observation 
that JLPZ gauge fixing can be "undone" by multiplying Eq. (H) (for = 0) by a con- 
stant in the form / Vnd, and transforming A^ A^, removing S^i from the partition 
function. Wave-function renormalizations are not necessary for physical quantities. 

In other words, renormalized perturbation theory for physical quantities is the 
same in both versions for a suitable definition of the renormalized coupling constant, 
but the relation between the renormalized and bare coupling constants is, in general, 
different in the FP and JLPZ versions of the theory. 

4-. In the case of noninvariant correlation functions, it is well known that, in 
the FP version of the theory, only a universal coupling-constant renormalization and 



5 



multiplicative wave-function renormalizations are necessary as a consequence of BRST 
invariance. The situation is less clear in the JLPZ version of the theory. Some 
progress can be made however, by observing that a set of local and renormalizable 
Feynman rules can be constructed for the JLPZ partition function, Eq. (0). Start 



with diagonalizing the quadratic form in 6 and rj in Eq. (|T5D by a shift 6 = 6' — rj. 



After this shift, Eq. can be written as the sum of two parts, 

S^A'^) = S^iA', + Si,{A', V, e') , (22) 



with Sl:^{A\ 1], 6' = 0) = 0. To order g/M we obtain 



upon using d^A' ^ = 0. Next, we define the following action (dropping the primes on 



Afj, and 



Crf = Si^,{A + dv/M) + S^,{A + dv/M)-S'l,{A,r]) + S'^,{A,r],e) (24) 

+^d,9%D,{A)9r - gr'^'d.e^^v'' + O 

Here ^^0{g/M)" indicates all higher order terms in the expansion in g, including 
those coming from the shift Eq. (p^ , and they should be kept (up to the order in g 
of interest). This action will give rise to a correct set of Feynman rules if we add the 
rule that a factor —1 be applied for each connected 6'-sub diagram without external 6 
lines. This is reminiscent of the minus sign for ghost loops in the FP case: the integral 
over 9 in the denominator of Eq. (^) gives rise to an effective action Ses{A,ri), 

exp(5eff(A,r/)) = Jve exp {-S'^,{A,r^,9)) , (25) 

the vertices of which correspond precisely to these connected ^-subdiagrams. The 
additional minus sign corresponds to the fact that this integral appears in the de- 
nominator of the integrand in Eq. (^. If we work with a regulator in which it is 
important to keep the nonlinear terms in the measure, this can be taken into account 
by treating these nonlinear terms as part of the action. 

Finally, there is still the 5-function of Eq. ([T0|), through which the field t] was 
introduced. After the shift Eq. (|T^, this is just 5(9^ A^), which again can be repre- 
sented as in Eq. ([191) . This leads to a set of Feynman rules for a local theory, Eq. (p^), 
which is renormalizable by power counting for any value of ^, with the original JLPZ 
version (Eq. (^)) corresponding to the limit ^ — >■ (at fixed cutoff). 

It is instructive to compare our local Feynman rules with the non-local Feynman 



rules of Ref. |T^. Our Feynman rules are defined by Eq. (0), with the minus-sign rule 
for connected ^-subdiagrams, and the 5-function constraint on A^, enforced by taking 
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^ — »• in Eq. (p^ . Note that the longitudinal part of the gauge field is represented 
by rj, through Eq. (|TUp. The Feynman rules of Ref. were derived by integrating 
out the field 9, or equivalently h, in Eq. (^), and without introducing the field rj to 
represent the longitudinal part of the gauge field. This leads to non-local vertices, 
where the non-locality arises from the integration over 6, as well as from the non-local 
relation between rj and the longitudinal part of A^. 

It is straightforward to obtain these non-local Feynman rules from our Feynman 



rules by (perturbatively) integrating out 6 in Eq. (^41) , and by converting the Feynman 
rules for rj to Feynman rules for the longitudinal part A^^ of A^. This is obvious, 
because both our Feynman rules and those of Ref. |jl2| are derived from the same 
partition function, Eq. (^). In momentum space, the relation between rj and the 
longitudinal part Al^ of A^ is 



Vip) = -tM^AL.ip) , 
and from {ri{p)r]{q)) = S{p + q)/p^ we thus find 



(26) 



Pf^P 



(27) 



in accordance with the first equation in Eq. (pTD, as well as with Ref. |]T2|. The 
simplest non-local vertex of Ref. [|l^ corresponds to the drjA t] three-point vertex of 
Eq. (p^, which reads in momentum space, with k the momentum of and p, q the 
1] momenta, 

'%p,-q,)- (2^ 



Using Eq. (E^), this translates into a non-local three-point vertex for the gauge fields: 
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q\k. 



KPf^ - qt^)^^ + KqK - ^k)^7^ + ^(^A - Pa) 



p2q2 



(29) 



where we symmetrized in the three gauge-field lines 
three-point gauge-field vertex of Ref, 
rules of Ref. p 
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This is exactly the non-local 
The other non-local vertices of the Feynman 



can be obtained in a similar way. 
It follows that, for correlation functions involving only gauge fields on the external 
lines, one obtains the same result from either our Feynman rules or those of Ref. ||12|| . 
This is true in particular for the one-loop vacuum polarization, which was calculated 
explicitly in Ref. [|r^. The result found there agrees with the vacuum polarization in 
Landau gauge, i.e. as calculated from Eq. (|20|). 

We conclude this section by noting that other field redefinitions of the field 6 can 
be used, for instance 



exp \ ij^e^ 



exp ( ij^e' ) exp 



or 
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exp ( i^6' I = exp (i^6'' I (l — ^vr'?) • 



(30) 
(31) 
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These examples differ only by terms of order g/M from the one employed above, and 
thus lead to a different specific form of the 0{g/M) terms in Eq. (0). 



5. We argued that the gauge-fixing procedure proposed in Refs. [T^, with 
the choice of Eq. @ for Sni{A), is perturbatively equivalent to the standard gauge- 
fixing procedure with Fadeev-Popov ghosts. In our derivation, we chose M to be of 
the order of the cutoff, thus extending earlier arguments in which the (formal) limit 
M — ^ oo was considered. 

"Equivalent" here means that perturbatively calculated relations between physi- 
cal quantities will be the same in both versions of the theory. Since, in addition, the 
JLPZ method leads to a weighting of the integration over orbits which takes Gribov 
copies correctly into account, this method may be the preferred one for nonperturba- 
tive calculations in gauge-fixed Yang-Mills theories. 

We also derived local Feynman rules for the JLPZ version of the theory, Eq. (Q), 
from which it can be seen that the theory is renormalizable by power counting. The 
choice of M at the order of the cutoff is a key ingredient here. By construction, 
correlation functions of gauge-invariant operators are the same when calculated per- 
turbatively from either Eq. (^) or Eq. (0), after the field redefinition of Eq. (|l^) is 
taken into account. In order to renormalize correlation functions of gauge noninvari- 
ant operators, counterterms may have to be added to Eq. (|2^), in addition to those 
needed to renormalize gauge-invariant quantities. The locality guarantees that all 
counterterms necessary for renormalization are local, and the renormalizability guar- 
antees that only a finite number, all with mass dimension less than or equal to four, 
will be needed. It is also clear that this can be done in such a way that the invariance 
of gauge-invariant correlation functions is maintained, because of the gauge invari- 
ance of the original formulation of Eq. (^. Hence, we believe that no problems will 
be encountered in carrying out this program order by order in perturbation theory, 
for the theory of Eq. (0). But it is not clear how this would then "translate back" 
to a JLPZ- like formulation as in Eq. (Q). What is lacking is a tool similar to BRST 
symmetry in the FP version of the theory, which could be used to further control the 
form of the counterterms. It would be interesting and useful if such a mechanism 
could be found. 

We end with a comment on our use of the specific form of the action 5*^(^4) in 
Eq. (^. This choice is the "most" (and only) relevant local operator in the sense 
of the renormalization group. Our analysis does not work when 5*111 (A) is chosen to 
be a marginal operator. We expect that they will work if we would add a marginal 
operator to Sni{A) of Eq. (|^). For instance, a term of the form cX)/^ tr(74^), with c 
a constant of order g"^ (which is natural on the lattice), can be removed by a field 
redefinition of the form — > + (c/2)A^/M^. Since the nonlinear term of this 
field redefinition is of order this will just remove the term X]^tr(y4^), without 

introducing any other marginal terms. 
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